Atypical Fractional Quantum Hall Effect in Graphene at Filling Factor z/g = 1/3 
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We study, with the help of exact diagonahzation calculations, a four-component trial wave function 
that may be relevant for the recently observed graphene fractional quantum Hall state at a filling 
factor UG = 1/3. Although it is adiabatically connected to a 1/3 Laughlin state in the upper spin 
branch, with SU(2) valley-isospin ferromagnetic ordering and a completely filled lower spin branch, 
it reveals physical properties beyond such a state that is the natural ground state for a large Zeeman 
effect. Most saliently, it possesses at experimentally relevant values of the Zeeman gap low-energy 
spin-flip excitations that may be unveiled in inelastic light-scattering experiments. 
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The recent observation of the fractional quantum Hall 
effect (FQHE) in graphene \m,\M, has proven the relevance 
of Coulomb interactions in this novel two-dimensional 
(2D) electron system, in agreement with theoretical ex- 
pectations [3|-|8[. The most pronounced state is the one 
found when the ratio vq = n^ijnB between the electronic 
density ng; and that of the flux quanta tib = eB/h is 
1/(3 = 1/3. Although this state is reminiscent, at first 
sight, of the prominent 1/3 state observed in semicon- 
ductor heterostructures [9|, which is described to great 
accuracy by the Laughlin state [lO|, several questions 
arise when taking fully into account the four-component 
strucure of graphene, due to its four-fold spin-valley de- 
generacy. Whereas first numerical approaches [5| con- 
sidered the physical spin to be frozen by the Zeeman 
effect and concentrated on the valley-isospin degree of 
freedom in a two-component system, a four-component 
approach [8| seems to be more appropriate in view of the 
rather small energy scale associated with the Zeeman ef- 
fect Az, when compared to the leading energy scale of 
the Coulomb interaction, e^ / elB at the magnetic length 
Ib = y/hc/eB. Indeed for a 5-factor of 2 ll| , one obtains 
Az/ie'^/elB) ~ 0.002^B[T] x e, where e is the relative 
dielectric constant. 

A further complication arises in graphene, as compared 
to the 2D electron gas in semiconductor heterostructures, 
when one considers the definition of the filling factor 
I'G, which is proportional to the carrier density Uei- In 
graphene, the carrier density vanishes at the Dirac point, 
where the spectrum is particle-hole symmetric. In the 
presence of a magnetic field, a four-fold degenerate zero- 
energy Landau level (LL) is formed that happens to be 
half-filled when Uei = and thus i^c = - the situa- 
tion at 1/(3 = is therefore more reminiscent of a fill- 
ing factor of ^/ = 2 in a usual four-component quantum 
Hall system [8[, and the observed FQHE at vq = 1/3 
corresponds to a situation where two of the four spin- 
valley subbranches are completely filled and a third one 
1/3-filled {u = 2 + 1/3). As a consequence the ob- 



served FQHE is not a simple Laughlin state with an 
SU(4)-spin-valley ferromagnetic ordering, which would 
arise at vq — ~2 -I- 1/3 (or by particle-hole symme- 
try, at i/(3 = 2 — 1/3) [y, |8|. A natural candidate for 
large values of the Zeeman gap would then be a valley- 
SU(2)-ferromagnetic Laughlin state ^V ^ in the spin- 

^ branch of the zero-energy LL similar to the usual 1/3 
physics. In this scenario, both states K and K' are com- 
pletely filled in the spin-f branch. The small relative 
value of the Zeeman gap, however, casts doubts on such 
a scenario of complete spin polarization induced by an 
external field, without considering a cooperative effect 
mediated by the Coulomb interaction. 

Here, we analyse the system in the zero-energy LL 
with the help of exact-diagonalization (ED) calculations 
for relativistic electrons in the spherical geometry with 
SU(4)symmetry that interact via the Coulomb interac- 
tion [J, la y. We show that already for a very small 
Zeeman effect, one may obtain a FQHE at i/q — 1/3 in 
graphene. This state may be described in terms of a four- 
component Halperin wave function ^2-1-1/3 '^hich is adi- 
abatically connected to the valley-SU(2)-ferromagnetic 
state in the upper spin branch. The latter is the natural 
ground state for a large Zeeman splitting. Most saliently, 
in spite of this adiabatic connection, the low-energy ex- 
citations in an intermediate range of the Zeeman split- 
ting are different from those of the simple SU(2) Laugh- 
lin state. In addition to the charge-density-wave mode 
with its characteristic magneto-roton minimum and the 
valley-isospin wave, which is the Goldstone mode asso- 
ciated with the spontaneous valley-isospin breaking in 
the spin- 4, branch, we find a low-energy spin-flip mode 
with a gap that depends linearly on the Zeeman cou- 
pling. These modes may be experimentally accessible 
in inelastic light-scattering measurements that have re- 
vealed similar modes in conventional quantum Hall sys- 
tems in GaAs heterostructures 
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That electrons 
in graphene reside at the sample surface makes this novel 
2D electron system even better adapted to optical mea- 



surements than the latter. 

In order to describe the FQHE state at uq — 1/3, 
which corresponds to a filhng factor oi v ~ 2 + 1/3 when 
counted from the bottom of the central n = LL, we 
investigate the trial four-component wave function 
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where z^ denotes the complex coordinate of the j'-th 
electron in the spin- valley subbranch a, ^ {a ^t or i a-nd 
^ = K 01 K'). We have omitted an ubiqitous Gaussian 
factor in the expression. Notice that, in the absence of 
a symmetry-breaking field, the wave function ([T|) is not 
a good trial state because the Coulomb interaction po- 
tential respects the SU(4)-spin- valley symmetry j4l Il4|. 
whereas the wave function ([Ij is not an eigenstatc of the 
SU(4)-Casimir operators. 

This is indeed corroborated by our ED calculations 
for A^ = 17 particles with Nb = 6 flux quanta thread- 
ing the sphere, the relation between N and Nb being 
Nb = (3/7) A^ - 9/7 for the state ^ [H, which yields 
the required filling v = 7/3 in the thermodynamic limit. 
The ground state is then found in spin sectors differ- 
ent from that, 25*2 = 11, expected for the state dD). A 
simple manner to stabilize the state ((!]) is to use appro- 
priate pseudo-potentials [l6[ that break the SU(4)-spin- 
valley symmetry in the interaction potential. However, 
surprisingly, this trial state becomes the ground state 
also when the SU(4) symmetry is broken by an exter- 
nal field - e.g. a very small value of the Zeeman effect 
turns out to be sufficient, A^ — O.Ole^/eZs, which is a 
tiny fraction of the leading Coulomb interaction energy 
scale e^ /els ^ 625(-\/i?[T]/e) K. In a typical experimen- 
tal situation, the 1/3 state has been found at a magnetic 
field of roughly B ~ 9. ..12 T [ll, I4I, which corresponds 
to a ratio of l^zKe^/dB) ^ 0.006. ..0.008 x e if one uses 
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FIG. 1: (Color online) Energy spectrum for A'^ = 17 electrons 
at a filling factor vq = 1/3 {v = 2 + 1/3), as a function of 
A2, obtained from ED calculations of the Coulomb interaction on 
the sphere {Ng = 6) with implemented SU(4) symmetry. Above 
A^ — O.Ole^ /cIb, the ground state is found in the maximally spin- 
polarized sector {2Sz = 11, red diamonds). The inset shows a zoom 
on the region for small values of A^. 
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FIG. 2: Classification of the excitations of the 'I' 
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excitations of a one-component Laughlin state are found in the 
same spin- valley sector (C), whereas the Goldstone mode due to 
the broken SU(2) valley-isospin symmetry in the spin-J, branch is 
an insospin-wave mode (ISW). In addition to these conventional 
modes, the four-component state |[TJ possesses a spin-flip (SF) 
mode. 



man gap Az- Above the critical value A^, the ground 
state is found in the maximally-polarized spin sector 
that corresponds to the state ([T]), whereas the excited 
state with the lowest energy is in the same spin sector, 
2Sz = 11, only above a second value A| ~ 0.03e^/eZB. 
g ^ 2 [ll|. This value is slightly smaller than our the- For values of the Zeeman gap A^ < A^ < A|, the 

excited state with lowest energy is found in the spin sec- 
tor 2Sz = 9. Above A^, however, the energy cost of 
this spin-flip excitation (SF, see Fig. [2|) is larger than 
the lowest-lying excitation in the fully polarized sector 
2Sz = 11 (C in Fig. [2]). 

These results suggest that the state ([T]) may have phys- 
ical properties beyond the simple 1/3-Laughlin state in 
the spin- 4, branch, in the form of coherent spin-flip excita- 
tions in an intermediate range of Zeeman gaps. In order 
to test this scenario in more detail, we have investigated 
the two-component wave function 



oretical estimate if one considers the smallest possible 
value of the dielectric constant (e = 1 for free-standing 
graphene). However, virtual interband excitations lower 
the dielectric constant that becomes Cno ~ 4 for free- 



standing graphene in the large- wave- vector limit 17[ 1 a-lso 
in a strong magnetic field [l8|, and the precise value of 
the dielectric constant in graphene remains an open issue. 
Notice that it is even still under debate whether the Zee- 
man splitting is indeed the dominant SU(4)-symmctry- 
breaking perturbation or whether the valley splitting is 
more relevant. Our theory and the conclusions of this 
paper, however, remain valid also in the latter case if one 
interchanges the role of spin and valley isqspin and if one 
replaces Az by a "valley Zeeman effect" [l9j . 

The energy spectrum, which we have obtained in ED 
calculations, is shown in Fig. [T]as a function of the Zee- 
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which would be a candidate in a two-component quantum 



Hall system, such as a conventional 2D electron gas in a 
GaAs heterostructure, at a filling factor v = 1 + 1/3. It is 
insofar related to the four-component wave function ([T]) 
as it describes the same physical situation if the valley- 
isospin degree of freedom for spin-| electrons is neglected. 
The novel wave function ^ therefore does not reveal 
any valley-isospin-wave mode (ISW, see Fig. [2|) that is 
the Goldstone mode of the spontaneously broken valley- 
SU(2) symmetry in the spin-J, branch and that may even- 
tually become gapped if one takes into account a possible 
valley splitting. In contrast to its four-component ana- 
logue ((!]), the two-component wave function ([2]) allows for 
a more detailed study of different system sizes in ED with 
an implemented SU(2) symmetry. Indeed, our ED cal- 
culations with an implemented SU(4) symmetry allowed 
only for one single system size {N = 17 particles, only 
^i = 3 populate the upper spin branch), in which case 
the subspace with maximal spin polarization {2Sz = 11) 
is of dimension one such that the overlap with the wave 
function ([T|) is trivially 1. 

Figure El^a) shows the energy spectrum for N ~ 22 
particles and Nb = 15 flux quanta, in the different spin 
sectors, obtained by ED of the SU(2) Coulomb inter- 
action potential in the lowest LL. The spectrum is in 
qualitative agreement with that obtained for the four- 
component system at t/ = 2 -I- 1/3 (Fig. [T]) - because 
the wave function ([2]) is not an eigenstate of the SU(2)- 
symmetric Coulomb potential, it does not describe the 
ground state at Az = 0, where one obtains a three-fold 
degenerate state (with 2Sz = 0,±2), but in a compress- 
ible sector (L ^ 0). As for the four-component case, a 
small symmetry-breaking Zeeman gap A^ — O.Ole^/e/s 
suffices to stabilize a state with maximal spin polariza- 
tion {2Sz = 10 and N^ = 6), which has an overlap of 99% 
with the wave function ([2]) [23|, and the lowest- lying ex- 
cited state in an intemediate range of the Zeeman gap, 
A^ < Az < A|- ~ 0.08e^/elB, involves a spin flip as it 
is found in the spin sector 2Sz = 8. 

It has been argued that, for vanishing Zeeman split- 
ting, the state at i^ = 1 + 1/3 should be a spin-singlet 
composite-fermion (CF) state with reversed flux attach- 
ment l20|. Hund's rule, according to which the system 
chooses a maximally polarized spin inside each energy 
level, would then predict an unpolarized state because 
1/ = 2/3 corresponds to a completely filled lowest CE- 
LL [20| . but the same rule favors a completely polar- 
ized state if applied to the original electron coordinates. 
Our results indicate that already for a very small Zee- 
man splitting, a completely polarized state is favored 
that satisfies the electronic instead of the CF version 
of Hund's rule. Notice, however, that a direct numeri- 
cal comparison between both states, CF spin singlet and 
state ([2]), is problematic in the spherical geometry be- 
cause the spin-singlet state has a different flux-particle- 
number relation, Nb = (3/4)A^ — 1, than the polarized 
state ©, TVs = (3/4) TV - 3/2. We find for TV = 20 and 
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FIG. 3: (Color online) (a) Energy spectrum for N = 22 electrons 
at a filling factor u = 1 + 1/3, as a function of Az in the differ- 
ent spin sectors 2Sz , obtained from ED calculations of the Coulomb 
interaction on the sphere {Nb = 15) with implemented SU(2) sym- 
metry. The inset shows a zoom on the region for small values of 
Az- (b) Excitation spectrum at Az = 0.05e /eZs, as a function 
of the angular momentum L. The energy is measured with respect 
to the ground state (at L = 0). The spin-fiip mode is found in 
the 2Sz = 8 sector (blue squares) and scales linearly in B, whereas 
the 2Sz = 10 sector reveals the usual magneto-roton branch (red 
diamonds), which scales as vB with the _B-ficld. 



Nb — 14 (results not shown) that the ground state is 
indeed a singlet at low Zeeman splittings, but it is max- 
imally polarized above a value of Az/{e^ /cIb) '^ 0.03, 
which is on the same order of magnitude as A^. 

In order to gain further insight into the nature of 
the low-lying excitations, we have calculated the spec- 
trum [Fig. El^b)] at an intermediate value of the Zee- 
man gap, Az = O.Obe"^ /elB, where spin- flip excitations 
are expected to be relevant. The spectrum is now plot- 
ted as a function of the angular momentum in order to 
make apparent possible low-energy collective excitations 
of the incompressible state ([2]). Within the charge sec- 
tor with no change in the spin polarization, one observes 
in Fig. pfb) the usual magneto-roton branch (red dia- 
monds) 12 ll which arises from gapped density- wave exci- 
tations J22l | and which is a prominent feature of Laughlin- 
type physics. However, another mode is apparent in Fig. 
^h) that indicates the presence of collective excitations 
beyond the usual one-component Laughlin state and that 



is precisely a spin- flip excitation (blue squares). This 
mode, which is the lowest-energy excitation in the low- 
L regime, is well separated from the high-energy part 
of the excitation spectrum, such that it is likely to be a 
true collective mode. Notice that in the large-L limit, the 
magneto-roton branch has a lower energy, and one may 
thus conjecture that the activation gap, i.e. the energy 
to create a well-separated quasiparticle-quasihole pair at 
large values of L, does not involve a spin-flip excitation, 
but is governed by one-component Laughlin physics. 

The relevance of collective spin-flip excitations in an 
intermediate Zeeman-gap range may eventually be tested 
experimentally in inelastic light-scattering experiments 
that are capable of probing collective excitations at finite 



wave vectors 



12 



13| . Indeed, these experiments probe 
characteristic parts of the dispersion relation that show 
an enhanced density of states (such as at its minima and 
maxima). Because the spin-flip mode in Fig. [3l^b) is 
almost flat at low angular momenta L that correspond to 
small wave vectors, one may expect an enhanced peak in 
such inelastic light-scattering measurements, at energies 
in the Q.le^ /els range (roughly half of the energy of the 
magneto-roton minimum, for the particular choice A^ = 
O.OSe^/e^s). As one may see in Fig. |3l the spin- flip 
excitation scales linearly with the Zeeman gap, such that 
the associated peak is expected to scale linearly with the 
magnetic field as well, whereas that of the usual magneto- 
roton would scale as ^/B [see Fig. [3l^b)]. The observation 
of such a linear S-field dependence of the light-scattering 
peak would be clear evidence for the relevance of spin- 
flip, beyond the properties of the Laughlin liquid, of the 
VG = 1/3 state in graphene. 

In conclusion, we have shown, within ED calculations 
for a four- and a two-component system on the sphere, 
how a FQHE can arise in graphene at vq = 1/3 even at 
very small values of a spin- valley symmetry-breaking Zee- 
man field. Although the leading energy scale is given by 
the SU(4)-invariant Coulomb interaction, a small Zeeman 
gap /S.z/[e^ /els) ~ 0.01 is sufficient to fully polarize the 
electronic spin and thus to stabilize the state ^ which 
we have identified as being responsible for the observed 
graphene FQHE [l|, l2| . In spite of its reminiscence with 
the Laughlin state, novel collective excitations that are 
inherent to the four-component character of graphene de- 
termine the low-energy spectrum at intermediate values 
of the Zeeman gap, A^ < A^ < A|, that correspond to 
the experimental situation in which the FQHE has been 
observed. In order to gain further insight into the nature 
of these spin-fiip excitations, which may be visible in in- 
elastic light-scattering experiments, we have performed 
ED calculations in an analogous two-component quan- 
tum Hall system at a filling factor v = 1 + 1/3 that 
corresponds to a completely filled spin-f and a one-third 
filled spin- J, branch. The spin-flip excitation is well sep- 
arated from the high-energy part of the energy spectrum 
thus indicating its collective nature, in addition to the 



usual magneto-roton branch that determines the low- 
energy spectrum in the large angular-momentum regime. 
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